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WALL  PRESSURE  EVENTS  FROM  THE  DIRECT  NUMERICAL  SIMULATION 
OF  A  TURBULENT  WALL-BOUNDED  FLOW 


INTRODUCTION 

Although  wall  pressure  events  beneath  a  turbulent  boundary  layer  have  been  expaimentally 
detected  with  inessure  transducers  mounted  in  the  wall,  the  number  of  spatial  locadtxis  at  which 
the  wall  pressure  can  be  measured  is  limited.  The  wall  pressure  frequency  spectrum  or  wall 
pressure  events,  such  as  pressure  peaks,  can,  however,  be  measured  at  a  handful  of  locations  as  a 
fiincticm  of  time,  although  the  character  of  the  pressure  field  in  space  cannot  be  determined  from 
these  measurements  at  a  few  discrete  locations.  On  the  other  hand,  direct  numerical  simulation 
of  turbulent  wall-bounded  flow  o^rs  access  to  the  wall  pressure  field  across  an  entire  region  of 
the  flow,  permitting  the  calculation  of  the  wavenumber  spectrum  and  the  detection  of  wall 
pressure  events  in  space.  Unfortunately,  the  nature  of  direct  numerical  simulations  is  such  that 
the  Reynolds  number  is  artificially  low  and  the  duration  of  time  of  the  simulatitm  is  severely 
limited.  Nevertheless,  wall  pressure  data  from  direct  numerical  simulatim  of  a  turbulent  channel 
flow  provide  a  spatial  representation  of  a  turbulrat  field  that  is  otherwise  unattainable. 

Conditional  sampling  of  experimental  wail  pressure  data  has  been  used  to  characterize 
significant  pressure  "events'*  related  to  the  turbulent  velocity  field.  Typically,  the  wall  pressure 
and  the  velocity  are  measured  simultaneously  at  discrete  locations  in  the  flow,  and  events  are 
detected  with  a  variety  of  techniques.  Use  of  this  process  in  a  planar  boundary  layer  has  allowed 
positive  peaks  in  the  wall  pressure  to  be  cmrelated  witii  streamwise  acceleration  related  to  near¬ 
wall  shear  layers.  Since  such  shear  layers  are  closely  related  to  turbulence  generation 
mechanisms  near  the  wall,  these  measurements  indicate  that  positive  wall  pressure  peaks  are 
coupled  with  turbulence  production.  Apart  from  conditionally  sampled  events,  though,  positive 
wall  pressures  generally  are  associated  with  high  streamwise  velocity,  and  negative  pressures  are 
associated  with  low  streamwise  velocity In  a  boundary  layer  on  a  cylinder,  large-amplitude 
negative  pressure  peaks  have  been  associated  with  local  decreases  in  streamwise  velocity, 
indicating  a  bidirectional  relationship  or  coupling  between  large-amplitude  wall  pressure 
fluctuations  and  the  rate  of  change  of  streamwise  velocity  near  the  wall.^ 

Results  are  reported  here  for  the  detection  of  wall  pressure  events  in  the  National  Aeronautics 
and  Space  Administration/Ames  Flight  Center  (NASA-Ames)  direct  numerical  simulation  data 
base  of  turbulent  channel  flow.  The  goal  of  the  work  was  to  characterize  the  nature  of  pressure 
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peak  events  in  die  spatial  domain.  Unfortunately,  velocity  data  are  not  available  fm  this  data 
base,  and  so  the  velocity  field  and  the  inessure  at  the  wall  caiuuH  be  related  to  one  anmher  as  in 
the  studies  described  in  references  1  through  S.  In  spite  of  this  shcntcoming,  the  numerical  data 
base  {ffovides  infonnation  on  die  spatial  character  of  wall  pressure  events  that  cannot  be  obtained 
fitmn  experiments  with  transducers  at  a  limited  number  of  points.^*^  In  the  work  presented  here, 
the  development  of  wall  pressure  events  in  both  time  a/id  space  is  invesdgared.  The  results  help 
provide  an  understanding  of  the  turbulent  nature  of  the  wall  pressure  fields  in  wall-bounded 
turbulent  flow. 


NUMERICAL  DATA  BASE 

A  direct  numerical  simulation  of  turbulent  channel  flow,  described  in  detail  by  Kim  et  al.,^ 
was  used  to  generate  the  pressure  statistics  presented  in  this  reptMrt  This  data  base  is  the  same  as 
that  used  by  Kim'^  and  by  Choi  and  Moin^  for  analysis  of  the  structure  and  space-time 
characteristics  of  turbulent  wall  pressure  fluctuations.  The  computational  domain  for  the 
simulation  is  Lx  X  Ly  X  L2  of  4x5  X  25  X  (4xf3)5  in  the  streamwise  x,  wall-n(xrmal  y,  and 
spanwise  z  directitms,  where  5  is  the  half-channel  width.  The  number  of  grid  points  N  and  other 
details  of  the  computational  domain  and  flow  parameters  are  provided  in  table  1.  The  simulation 
was  carried  out  for  25,640  time  steps,  saving  tire  pressure  for  the  entire  plane  of  one  of  the  walls 
at  every  tendi  step  corresponding  to  At  »  3.75  x  10*3S/Ut.  The  data  base  was  stored  in 
MATLAB-compatible  format  on  a  Sun  workstation  after  translation  from  the  original  Cray- 
compatible  stOTage  format  The  results  described  here  are  based  on  the  first  1500  time  steps, 
corresponding  to  a  total  time  of  5.625  5^f 

Table  1.  Parameters  of  the  Direct  Simulation  Data  Base^ 


Res  =  Uci^=  3200 
Ree=Ucj0A^  =  280 
Ret*Ut5^^  =  180 

Streamwise: 

Lxy6*4x 

Nx=128 

Ax+=17.6  Lx+  =  225: 

Wall-normal: 

Ly/5  —  2 

Ny=129 

Ay+  =  0.05  (near  wall)  to 
Ay+  =  (centerline) 

Spanwise: 

Lz/5=4x/3 

Nz=128 

Az+  =  5.9  U+  =  755 

AtUt/5=  3.75x10-3 
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AlthcMJgh  this  numerical  data  base  provides  a  vast  set  of  wall  pressure  data,  there  are  several 
^wbacks.  Hrst,  the  momentum  thickness  Reynolds  number  for  the  simulaticm,  Ree  =  280,  is 
probably  physically  unrealizable,  although  the  low  Reynolds  number  is  necessary  because  of 
computational  limitarjAny.  Second,  periodic  boundary  conditicMis  require  that  the  fluid  stroctures 
entering  the  conqtutational  domain  be  identical  to  those  leaving  tlw  dmnain,  which  cannot  occur 
in  a  real  channel  flow.  Third,  the  extent  of  the  computational  domain  is  very  limited,  2253 
viscous  units  in  the  streamwise  direction  and  755  viscous  units  in  the  spanwise  directitm.  Fourth, 
tlw  wavenumber  spectrum  of  the  wall  pressure  has  a  suspicious  increase  in  energy  as  the 
wavenumber  approaches  zero.^  Finally,  the  p(»id<m  of  die  centerline  of  tlw  simulated  channel  is 
at  a  very  small  distance  from  the  wall  corresponding  to  5UxA^180,  which  suggests  that 
contributitMis  by  the  inner  region  of  the  flow  may  dominate  because  the  outer  flow  is  largely 
absent 


RESULTS 

A  contour  plot  of  the  wall  pressure  at  a  typical  instant  in  time  is  shown  in  figure  la.  Given 
the  spacing  for  the  contours  of  1  pnns«  it  is  evident  that  most  of  the  field  has  relatively  low  levels 
of  pressure  fluctuations.  Only  within  the  solid  omtours  does  the  wall  pressure  exceed  ±  1  pmu- 
The  positive  and  negative  contours  for  p  >  prms  are  shown  in  figures  lb  and  Ic,  respectively. 
Both  positive  and  negative  contours  are  each  fairiy  uniformly  distributed  in  space  and  seem  to 
have  no  preferential  elongation  or  sh^.  But  positive  and  negative  pressure  regions  seem  to 
often  be  located  near  one  another  (for  instance,  at  z/5  ~  3  and  x/5  ~  5). 

Relatively  strong  pressure  events  are  tracked  in  time  in  figure  2a.  In  this  figure,  events 
exceeding  ±3  pmu  are  tracked  for  a  time  period  of  0.375  S/Ux  and  plotted  every  0.0375  S/Ut-  A 
series  of  similarly  sized  contours  in  a  vertical  line  shows  the  motion  of  the  pressure  event  Ftn* 
instance,  at  z/8  ~  1,  the  line  of  contours  extending  frcnn  x/5  ~  9  to  6  shows  the  translation  of  a 
negative  pressure  event  moving  from  top  to  bottom  in  the  domain.  A  typical  strong  pressure 
event  is  evident  at  z/5  "  3  and  x/5  ~  2  to  6.  The  streamwise  extent  of  the  event  is  several  5,  and 
the  duration  of  the  event  is  at  least  as  long  as  the  time  covered  in  the  figure,  which  is  0.375  S/Ut. 
Positive  and  negative  events  are  shown  in  Hgures  2b  and  2c,  respectively,  where  it  is  evident  that 
such  events  often  occur  near  one  another  (fOT  instance,  at  2^  ~  3  and  x/5  ~  2  to  6).  On  the  other 
hand,  positive  and  negative  events  do  not  necessarily  occur  near  one  another  (for  instance,  a 
negative  event  at  z/5  ~  1  and  x/5  ~  6  to  9). 
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Figure  1.  Contours  of  Wall  Pressure  at  a  Typical  Instant  in  Time 
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Substantial  infonnation  regarding  the  character  of  significant  wall  pressure  events  can  be 
obtained  from  the  analysis  of  "typical"  ones.  To  acccanplish  this,  pressure  peaks  exceeding  a 
trigger  level,  kpnas*  were  ensemble  averaged  to  gemrate  a  typical  event.  A  similar  technique  is 
used  in  analyzing  the  wall  pressure  time  signature  obtained  from  experimental  wall  pressure 
measurements.  However,  the  numerical  simulation  data  base  permits  the  detection  of  pressure 
peak  events  in  both  space  and  time.  The  detectiem  begins  with  a  search  for  peaks  exceeding  the 
trigger  level  for  a  single  instant  After  all  the  peaks  have  been  found  for  that  time  instant,  the 
next  time  instant  is  searched  for  peaks  exceeding  the  threshold.  Of  course,  most  peaks  found  at 
the  first  time  instant  are  also  found  at  the  second  time  instant  Each  peak  is  compared  to  the 
peaks  found  at  the  previous  time  instant  to  determine  which  peak  is  larger.  Subsequent  time 
instants  are  searched  for  peaks  arul  checked  against  nearby  peaks  from  previous  time  instants. 
As  this  process  continues,  the  location  and  time  instant  crmrespotKling  to  the  maximum  for  each 
pressure  peak  exceeding  the  trigger  are  found.  After  locating  all  pressure  peak  events  in  the  data 
base,  the  pressure  field  surrounding  each  detected  event  is  recovered  from  the  data  base  and 
ensemble  averaged  with  other  events  to  create  the  pressure  field  surrounding  a  typical  event 

The  results  of  such  an  analysis  for  positive  pressure  events  are  shown  in  figure  3.  Pressure 
contours  of  a  positive  event  for  k  s  4  are  shown  in  figure  3a.  The  pressure  event  has  been 
centered  in  the  figure  at  the  origin.  Although  most  of  the  surrounding  pressure  field  has  been 
averaged  out  to  nearly  zero  pressure,  just  upstream  and  downstream  of  the  positive  pressure 
event,  negative  pressure  occurs.  The  negative  pressure  area  is  relatively  weak  compared  with  the 
positive  pressure  event  The  negative  pressure  area  does  not  extend  to  the  side  of  the  pressure 
peak. 

Slices  through  the  positive  pressure  event  in  the  streamwise  and  spanwise  directions  are 
shown  in  figures  3b  and  3c,  respectively,  for  five  trigger  levels.  Of  course,  the  height  of  the  peak 
varies  with  the  trigger  level.  The  negative  pressure  just  upstream  and  downstream  of  the 
pressure  peak  is  evident  in  figure  3b.  The  upstream  negative  pressure  area  is  dependent  upon  the 
trigger  level.  Stronger  positive  pressure  peaks  result  in  stronger  negative  pressures  upstream  of, 
or  behind,  the  peak.  The  downstream  negative  pressure  region  is  weaker  and  somewhat  less 
dependent  upon  trigger  level.  The  pressure  peak  itself  has  a  streamwise  and  spanwise  extent  of 
about  80  wall  units.  The  negative  pressure  ahead  of  the  pressure  peak  extends  80  to  ISO  wall 
units  downstream,  and  the  negative  pressure  behind  it  extends  about  50  wall  units  upstream. 
Similar  results  are  evident  in  the  ensemble-averaged  pressure  for  Variable  Interval  Space 
Averaging  (VISA)  events  detected  based  on  the  streamwise  velocity.^  The  authors  used  the 
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Contour  Plot  of  a  Typical  Event  Exceeding  4  pn„s 
(BaUd  on  194  Individual  Events) 

Dotted  curves  are  zero  pressure  contours.  Solid 
curves  are  positive  contours  spaced  at  1  prms- 
Dashed  curves  are  negative  contours  spaced  at  025 
Prms-  Positive  pressure  peak  is  centered  at  the 
origin.  Upstream  corresponds  to  negative  x,  and 
downstream  corresponds  to  positive  x. 
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Cut  Through  the  Contours  Along  a  Streamwise 
Line  Through  the  Pressure  Maximum 

Trigger  levels  of  3,  4,  5,  6.  and  7  prms  ort 
identifiable  by  the  height  of  the  pressure  peak 
(largest  trigger  levels  are  highest). 


Cut  Through  the  Contours  Along  a 
Spanwise  Line  Through  the  Pressure 

Maximum 


Figure  3.  Ensemble  Average  of  Conditionally  Sampled  Positive  Pressure  Peak  Events 
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same  data  base,  but  only  peripherally  considered  the  pressure  at  the  wall.  Nevertheless,  they 
showed  that  a  positive  peak  in  pressure  is  related  to  a  shear  layer  directly  above  the  location  of 
the  maximum  pressure.  The  shear  layer  is  similar  to  positive  events  detected  experimentally 
with  Variable  Interval  Time  Averaging  (VITA).  Johansson  et  al.^  also  show  contours  of  pressure 
at  y-f  s  IS,  corresponding  to  a  VISA-detected  shear  layer.  The  negative  contours  at  this  distance 
from  the  wall  differ  considerably  from  those  at  the  wall  shown  in  figure  3a.  In  addition  to 
regions  of  negative  pressure  ahead  of  and  behind  the  pressure  event,  the  negative  contours  extend 
along  either  side  of  the  pressure  peak. 

A  similar  analysis  of  negative  wall  pressure  events  is  shown  in  rigure  4.  The  negative 
pressure  peak,  centered  at  the  origin  in  figure  4a,  lies  between  weak  positive  pressure  regions 
upstream  and  downstream  of  the  event  The  strength  of  the  positive  pressure  area  is  dependent 
upon  the  threshold,  as  indicated  in  figure  4b.  Higher  thresholds  result  in  a  higher  positive 
pressure,  especially  ahead  of  the  event  The  negative  pressure  events  are  not  as  broad  as  positive 
events  (see  rigure  4c). 

The  time  history  of  ensemble-averaged  positive  wall  pressure  events  is  shown  in  figure  Sa  for 
a  threshold  of  k  =  44.  Each  curve  is  the  wall  pressure  at  an  instant  in  time  along  a  streamwise 
line  aligned  with  the  spanwise  center  of  the  pressure  maximum.  No  wall  pressure  peak  is  evident 
in  the  lowest  curve  for  a  time  of  -0.225  5/Ut,  but  by  the  next  time  instant,  two  small  bumps 
appear  at  x***  =  -500.  As  time  progresses,  the  upstream  bump  grows  larger  than  the  downstream 
bump  (at  -0.15  5/Ut  and  -0.1125  S/Ux).  Then  the  downstream  bump  grows  much  larger  (at 
-0.075  SAJx  *tnd  -0.0375  8/Ux),  eventually  becoming  the  pressure  peak.  The  decay  of  the 
pressure  peak  shows  some  hint  of  the  double  pressure  peak  again  (at  0.1125  S/Ux  and 
0.1875  S/Ux).  The  pressure  peak  broadens  as  it  decays.  Similar  results  occur  for  trigger  levels  of 
k  =  3,  5, 6,  and  7.  Since  positive  pressure  peaks  have  been  associated  with  shear  layers  that  are 
related  to  the  bursting  process,  ^*4  it  is  reasonable  to  assume  that  the  double  bump  leading  to  the 
pressure  peak  is  also  related  to  this  process.  Howevor,  it  is  difficult  to  relate  the  double  bump  to 
common  conceptual  models  of  the  process  leading  to  a  burst,  such  as  a  haiipin  vortex. 

The  time  history  for  a  negative  event  is  shown  in  figure  5b.  In  this  case,  the  negative  bump 
leading  to  the  peak  is  preceded  by  a  slight  positive  pressure  (-0. 1125  ^x  through  -0.0375  S/Uxh 
This  leading  positive  pressure  becomes  the  relatively  strong  positive  pressure  ahead  of  the 
negative  pressure  peak.  After  the  negative  peak,  the  pressure  peak  broadens  and  decays. 
Negative  pressure  peaks  have  been  related  to  decelerations  of  the  streamwise  velocity.^  But  it  is 
unclear  if  negative  pressure  peaks  can  be  related  to  any  coherent  turbulence  generation  process  in 
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Contour  Plot  of  a  Typical  Event  Exceeding  -4  pnns 
(Based  on  230  Individual  Events) 

Dotted  curves  are  zero  pressure  contours. 
Solid  curves  are  positive  contours  spaced  at 
025  prms-  Dashed  curves  are  negative 
contours  spaced  at  1  prms-  Negative  pressure 
peak  is  centered  at  the  origin.  Upstream 
corresponds  to  negative  x,  ai^  downstream 
corresponds  to  positive  x. 


(a) 


Cut  Through  the  Contours  Along  a  Streamwise 
Line  Through  the  Pressure  Minimum 

Trigger  levels  of -3,  -4,  -5,  -6,  and  -7 prm  <tre 
identifiable  by  the  height  of  the  pressure  peak 
(largest  trigger  levels  are  Hghest). 


Cut  Through  the  Contours  Along  a 
Spanwise  Line  Through  the 
Pressure  Minimum. 


Figure  4.  Ensemble  Average  of  Conditionally  Sampled  Negative  Pressure  Peak  Events 
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ihgctf  *4 


Each  curve  is  a  cut  along  a  streamwise  line  through  the  pressure  peak.  Time  progresses 
from  the  bottom  curve  to  the  U)p  curve.  The  time  step  between  curves  is  0J0375  SlUr  (curves 
correspond  to  -0225.  -0.1875,  -0.15.  -0.1125,  -0.075.  -0.0375. 0, 0.0375, 0.075, 0.1125, 0.15, 
0.1875, 0225  from  bottom  to  top).  The  time  for  which  the  pressure  peak  occurs  is  the  seventh 

curve  from  the  bottom.  Flow  is  from  l^  to  right. 


H^uie  5.  Ensemble  Average  of  the  Time  Development  of  ConditionaUy  Sampled  Pressure  Peak 

Events  Exceeding  ±4  pnns 
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the  boundaiy  layer  in  the  same  way  that  positive  pressure  peaks  have  been  related  to  the  bursting 
process. 


CONCLUSIONS 

By  following  wall  pressure  events  in  both  space  and  time,  certain  characteristics  of  the 
evolution  of  these  events  have  been  found.  Positive  and  negative  events  often  appear  in  the  same 
localized  region  of  the  flow.  Ensemble  averages  of  conditionally  sampled  pressure  peak  events 
show  that  a  negative  pressure  appears  ahead  of  and  behind  a  positive  pressure  peak,  although  not 
at  the  sides  of  the  pressure  peak.  Similarly,  a  positive  pressure  appears  ahead  of  and  behind  a 
negative  pressure  peak,  but  not  at  its  sides.  Positive  pressure  peaks  have  been  associated  with  the 
shear  layer  related  to  bursting.  Based  on  the  results  of  Johansson  et  al.,^  the  streamwise  velocity 
is  negative  and  the  wall-normal  velocity  is  positive  above  the  negative  pressure  that  is  ahead  of 
the  pressure  peak,  and  vice  versa  for  the  negative  pressure  behind  the  pressure  peak.  As  a  result, 
the  small  negative  pressures  ahead  of  and  behind  the  positive  pressure  peak  are  related  to  high 
levels  of  Reynolds-stress-associated  production.^ 

The  relation  between  negative  pressure  peaks  and  flow  structures  has  not  been  discussed  to 
any  extent  in  the  literature.  Nevertheless,  negative  pressure  peaks  contribute  as  much  to  the  wall 
pressure  spectrum  as  do  positive  pressure  peaks,  are  as  energetic  as  positive  pressure  peaks,  and 
occur  as  firequently  as  positive  pressure  peaks.  Consequently,  a  full  understanding  of  the  wall 
pressure  requires  further  knowledge  of  how  flow  structures  are  related  to  negative  pressure 
peaks. 

The  sequence  of  events  leading  to  positive  pressure  peaks  includes  the  appearance  of  a 
double  positive  pressure  hump  before  the  maximum  pressure.  Eventually,  the  downstream  hump 
becomes  the  dontinant  pressure  peak,  which  then  decays  as  it  moves  downstream.  Without 
further  analysis,  it  is  unclear  how  the  double  hump  leading  to  the  pressure  maximum  is  related  to 
commonly  hypothesized  structures  and  events  in  the  flow.  The  pressure  signature  leading  to  a 
negative  pressure  peak  consists  of  a  negative  hump  behind  a  positive  hump.  The  negative  hump 
grows  substantially,  becomes  the  pressure  peak,  and  then  decays.  Again  it  is  difficult  to  relate 
the  development  of  the  flow  structure  to  a  flow  structure  in  the  boundaiy  layer. 
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